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Charged multiplicities in nucleus-nucleus collisions are calculated in the Dual 
Parton Model taking into account shadowing corrections. Its dependence on the number 
of collisions and participants is analyzed and found in agreement with experiment at 
SPS and RHIC energies. Using these results, we compute the J/ip suppression at SPS 
as a function of the transverse energy and of the energy of the zero degree calorimeter. 
Predictions for RHIC are presented. 



CHARGED MULTPLICITIES IN THE DUAL PARTON MODEL 



In the Dual Parton Model (DPM) the charged multiplicity per unit rapidity in a symmetric 
collision is given by ^ 
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Here P and T stand for the projectile and target nuclei, n{b) is the average number of binary collisions 
and nA{b) the average number of participants of nucleus A. These quantities can be computed in a 
Glauber model, k is the average number of inelastic collisions in pp and = kn{b)/nA{b) is the total 
average number of collisions suffered by each nucleon. The first term in (^) is the plateau height in a pp 
collision, resulting from the superposition of 2k strings, multiplied by ua- Since in DPM there are two 
strings per inelastic collision, the second term, consisting of strings stretched between sea quarks and 
antiquarks, makes up a total number of strings equal to 2kn. 

The charged multiplicity produced by a single string is obtained by a convolution of momentum 
distribution function and fragmentation functions (eqs. (3.1) to (3.4) of [^). 

Shadowing corrections in Gribov theory are universal ||^, i.e. they apply both to soft and hard 
processes. The reduction of the multiplicity resulting from shadowing corrections has been computed in 
^. These corrections are negligeable at SPS energies but at RHIC energies they reduce the multiplicity 
by 40 to 50%. 



We present the results ^ obtained at two different energies : ^/s = 17.3 and 130 GeV. Tlie 
corresponding non-diffractive cross-sections are (Jnd = 26 and 33 nib, respectively. We take k = 1.4 
and 2.0 corresponding to dN^J^ /dy = 1.56 and 2.72. The result in absence of shadowing at ^/s — 17.3 
is shown in Fig 1 (a). We obtain a mild increase of the multiplicity per participant consistent with 
the results of the WA98 Collaboration This increase gets stronger with increasing energies. As we 
pointed out before, shadowing corrections are negligeable at SPS energies but their effect is large at RHIC. 
Unfortunately, shadowing corrections have a rather large uncertainty at RHIC energies. Two alternative 
calculations [Q of shadowing lead to the results at ^/s = 130 GeV shown by the solid lines in Fig 1 (b). 
Clearly, with the larger values of the shadowing corrections we obtain a quantitative agreement with the 
PHENIX data [i. 



a) 



FIG. 1. a) The values of dN"'' /drj/upart versus npart for PhPb collisions at ^/s — 17.3 GeV in the range 
—0.5 < rjcm < 0.5 computed from eq. (§), compared with the WA98 data §. b) The values of dN"'' /dr]/(0.5npart) 
for Au-Au collisions at = 130 GeV in the range —0.35 < r/cm < 0.35 computed from eq. (0) taking into account 
shadowing corrections (see main text). They are compared to the data M. 



Note that our calculations refer to dN/dy while the first RHIC measurements refer to dN/d-q. 
The latter is, of course, smaller at mid rapidities. This difference is negligibly small as SPS where the 
laboratory pseudo-rapidity variable is used. However, at ^/s — 130 GeV where rjcm is used instead, their 
ratio can be as large as 1.3 

The predictions at yfs = 200 GeV were also given in ||l|. The predicted increase between 130 and 
200 GeV is 13 % in quantitative agreement with the measurement by the Brahms Collaboration Q. 



J/ij SUPPRESSION VS Et AND Ezdc AT SPS 

In the comovers approach the J/-0 survival probability is the product of two factors Sabsib, s) ■ 
Scoib, s). The first factor represents the suppression due to nuclear absorption of the cc pair. Its expres- 
sion, given by the probabilistic Glauber model, is well known. It contains a parameter, the absorptive 
cross-section dabs- The second factor Sco{b,s) represents the suppression resulting from the interaction 
with comovers. They are given by 



aabsf, X _ [1 - expi^ATAjs) (TgbM^ " exp(-5 Tsib-s) agbs)] , . 
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In Ny'^^{b, s) is the density of charged comovers (positives and negatives) in the rapidity region of 
the dimuon trigger and Nf — 1.15 fm~^ ||] |^ is the corresponding density in pp. The factor 3/2 in (^) 
takes care of the neutrals. In the numerical calculations we use aabs — 4.5 mb and CTco — 1 nib jlO| [Tl] ]. 
We compute the density of comovers in the framework of the DPM as we have explained before. 

This approach allows to compute the impact parameter of the J/^ event sample. Experimental 
results of this quantitie are plotted as a function of observable quantities such as Et or Ezdc- Using 
the proportionality between Et and multiplicity, we have 

ET{b)^^qN^:jb) . (4) 

The multiplicity of comovers (6) is determined using Eq. (|^) in the rapidity region of the NA50 
calorimeter (1.1 < yiab < 2.3). The factor 1/2 is introduced because Ny°^^{b) is the charged multiplicity 
whereas Et refers to neutrals. Thus the coefhcient q is close to the average energy per participant and 
its value can be determined from the position of the "knee" of the Et distribution of the MB events 
measured by the NA50 Collaboration. We obtain q = 0.62 GeV [|l2|. 
The energy of the zero degree calorimeter is defined as 

EzDcib) ^[A- nA{b)]E,r, + anA{b)E,^ . (5) 

Here A — n^(fe) is the number of spectator nucleons of A and Ein = 158 GeV is the beam energy. The 
last term represents the small fraction of wounded nucleons and/or fast secondaries that hit the ZD 
Calorimeter. The value of a can be precisely determined from the position of the MB event sample 
measured by NA50. We obtain a = 0.076. Eqs. (||) and (H) also lead to a correlation between 
(average values of) Et and Ezdc- This correlation is close to a straight line |jl^ and gives a good 
description of the experimental one. 

To explain the experimental data beyond the knee of the Et distribution we introduce two effects: 

1. Comovers fluctuations [|l^: We introduce the fluctuation in the density of comovers by 
replacing A^™^ in Eq. (|) by N^f^Tib, s)F{b) where F{b) = ET/ET{b). Here Et is the measured value 
of the trasverse energy and ETib) is its average value given by Eq. (^) - which does not contain the 
fluctuations. 

2. Et loss jl^: In the J/ij: event sample, Et ^ 3 GeV is taken by the J/^ trigger and, thus, the 
transverse energy deposited in the calorimeter by the other hadron species will be slightly smaller than 
the corresponding one in the MB event sample. 

The results of our model for the ratio J /-ip over DY , versus Et, in PbPb collisions at ^/s = 158 



GeV, are shown in Fig. 2a and compared with NA50 data |14 1^ - both for the true J /ip over DY ratio, 
and for the one obtained with the MB analysis. The results of the model for the true ratio are given by 
the dotted line (without fluctuations) and the dashed line (with fluctuations). In both cases the Et loss 
mechanism is not taken into account. However, our results for the true ratio J/t/i over DY do not change, 
since the effect due to the Et loss cancels in this ratio. We see that our results are in good agreement 



with the NA50 which do not extend beyond the knee. The other data in Fig. 2a are obtained with the 
MB analysis, and have to be compared with the dashed-dotted and solid curves (obtained taking into 
account the Et loss). In this case the agreement with the NA50 data is substantially improved. 

In Fig. 2b the results for the ratio J/tp over DY versus Ezdc are shown. These curves are 
obtained from the corresponding ones versus Et applying the Et — Ez dc correlation ||l3| . Comparing 
the data with the model predictions, we see that a better description of the central data (small Ezdc) is 
obtained when the Et fluctuations are taken into account. This was to be expected since the fluctuations 
in Et and Ezdc are related to each other via the Et — Ezdc correlation. It is important to note that 
the effect of the Et loss is not present in this case. Indeed Ezdc measures the energy of spectators and 
it is not affected by the dimuon trigger. No disagreement between the data and the model predictions is 
observed in this case for very central events. 
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FIG. 2. a) Ratio J/i) over DY versus Et in PbPb collisions at 158 GeV Q. The data are from ref. [0. The 
data labeled with DY are for the true J/ip over DY and they should be compared with the dotted and dashed 
lines obtained, respectively, without and with Et fluctuations. The data labeled Min. Bias should be compared 
with the dashed-dotted and solid lines, obtained with the Et loss, b) Ratio J/ip over DY versus Ezdc in PbPh 
collisions at 158 GeV per nucleon [^. The dashed (dotted) curve is obtained from the dashed (dotted) one versus 
Et applying the Et — Ezdc correlation and keeping the normalization unchanged. The data are from |l^ jl5|. 
The NA50 nuclear absorption curve is also shown. 



J/tp SUPPRESSION AT RHIC 

In order to compute the J /ip suppression at RHIC in the comovers approach we need to know not 
only the effects of nuclear absorption and comovers interaction, but also the effect of shadowing - which 
becomes important at RHIC energies. Therefore, our predictions are only valid so far as the shadowing 
corrections cancel in the J/V' over DY ratio - which is not necessarily the case [|l6| . 

For the comovers survival probability, Eq. (||), we keep the value CTco = 1 mb as discussed above. 
The hadronic multiplicity at RHIC energies have been successfully evaluated in DPM Therefore, the 
prediction at RHIC for this survival probability is rather safe. The situation is quite different in what 
concerns the survival probability due to nuclear absorption. It is widely recognized [pj|-pO| that, at high 



energy, when the coherence length becomes larger than the nuclear size, the probabilistic expression (g) 
is no longeer valid. It has been shown in that, at asymptotic energies, Eq. (||) is replaced by 

Sabs{b,s) = exp (^-i a A Ta{s)^ exp (^-i a B Tgib-s)^ (6) 

where a is the total cc-N cross-section. lia~ <7abs, the asymptotic result is not very different from the one 
obtained with the low energy formula - since Eqs. and (^) differ only in the second correction term. 
This situation is expected if the cc pair is produced in a colorless state interacting as a dipole. However, 
the possibility has been advocated Q , that the cc pair is produced in a color state accompanied by light 
quarks - in order to make the system colorless. In this case the system interacts with a comparatively 
large cross-section a ^ 15-^-20 mb 18 2^. With such a cross-section it is possible to explain the large 
suppression of the J/ijj observed at ~ 1, without initial state energy loss of the projectile partons in 
the nucleus. 
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FIG. 3. Ratio J/^p over DY versus charged multiplicity in Pb-Pb collisions at ^/s = 200 AGeV in the range 
—0.5 < y* < 0.5. The dotted line is obtained with the low energy nuclear absorption Eq. (^) and Gabs = 4.5 mb. 
The full line is obtained with the asymptotic expression Eq. and cr = 15 mb. 

The results for the ratio J/-0 over DY at mid-rapidities are presented in Fig. 3. The dashed line 
is obtained with Eq. (j^) and aabs — 4.5 mb. The solid line is obtained with Eq. (^) and a = 15 mb. As 
we see, the difference between the two predictions is very large. 

The proposed way of measuring the J /tp suppression at RHIC is via the ratio J /i}) over T. For 
the latter, the situation is even more complicated. Due to its larger mass, the correlation length is smaller 
and Eq. (||) is not valid even at mid-rapidities. In this case a finite energy formula [0110 ^as to be 
used - which interpolates between the low energy, Eq. (||), and the asymptotic (Eq. (^) limits. Actually, 
even for J/i/', Eq. (||) is not exact at mid-rapidities. We have estimated that using the finite energy 
formula with a = 15 mb the ratio J/ij] over DY for central events is 40 % higher than the solid line in 
Fig. 3. Therefore an accurate measurement in pA interactions will be necessary in order to clarify the 
theoretical situation. 
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